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Abstract. An algebraic derivation is presented which yields the exact solution of 
the mean first-passage and mean residence times of a one-dimensional asymmetric 
random walk for quenched disorder. Two models of disorder are analytically 
treated. Both, absorbing-absorbing and reflecting-absorbing boundaries are 
considered. Particularly, the interplay between asymmetry and disorder is studied. 
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1. Introduction 

The problems of the first-passage time (FPT) ^ [2j and the residence time (RT) [3] 
are very important issues in random walk theory. Moreover, several properties of 
diffusion and transport in disordered systems are based in this concepts. Here, we 
do not want to present a survey of the enormous literature in the field of mean first- 
passage time (MFPT) and mean residence time (MKT) of random walks, nevertheless, 
we wish to single out the works involved with analytical or exact results, mainly in 
one-dimensional disordered systems. Goldhirsch and Gefen |3] developed an analytical 
method for calculating MFPT for branched networks such as finite segments with 
dangling bonds and loops. The method is based on the generating function and was 
generalized for biased walks |S] ■ Extensions of the generating function method were 
done for analyzing the probability distribution function of FPT P| and the current 
autocorrelation function [7j. Later on, the generating function method was used for 
random one-dimensional chains [Hll], particularly for the Sinai problem ^U]. Explicit 
expressions for the MFPT, in terms of the basic jump probabilities for discrete time 
random walk with a reflecting boundary were obtained independently by Van den 
Broeck ^T], Le Doussal ^2): and Murthy and Kehr ^SI by different methods. It 
is interesting to remark that Gardiner JT1| had previously reported explicit MFPT 
formulae. An exact solution of the generating function for the first-passage probability 
was presented by Raykin ^2] using enumerative conbinatorics for summing up over 
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trajectories of the random walker jl5) . The distribution of escape probabihties was 
computed exactly by Sire J7] and the exact renormalization group analysis was 
performed by Le Dousal et at jl8| . In the last few years, one of the main applications 
of the exact expressions for MFPT and MRT in one-dimensional lattices was the study 
of exciton migration in treelike dendrimers (light harvesting antennae) HHI ■ 

A successful perturbative theory for survival statistics in disordered chains is the 
finite effective medium approximation (FEMA) ^Slj . Extensions of FEMA to biased 
media |221 and periodically forced boundary conditions were presented. A unified 
framework for the FPT and RT statistics in finite disordered chains with bias was also 
presented by the authors in reference j24| . where exact equations for the quantities 
averaged over disorder were obtained for both problems and its solutions up to first 
order in the bias parameter were constructed retaining the full dependence on the 
system's size and the initial condition. 

In this paper, we obtain explicit analytical expressions for the MFPT and MRT 
of random walks on a one-dimensional lattice for a quenched realization of disorder. 
Then, we consider two models for the disorder in the hopping transitions and average 
the expressions, in each case, on the realizations of disorder. The outline of the 
paper is as follow. The starting point of our formulation is given in section |21 In 
sections l2 . 1 1 and l2 . 21 we discuss the algebraic method that allows us obtain the exact 
dependence of the MFPT and the MRT on the set of transition probabilities. The 
biased nondisordered chain is treated in section|21 whereas two models of disorder with 
asymmetric transition probabilities are analyzed in sections 0] and The effects on 
MFPT and MRT of cutting the chain, at a reflecting site, are considered in sectionlHl 
Finally, in section 13 we briefly summarize the principal results of our study. 

2. Survival and residence probabilities 

We start considering a random walk on a discrete one-dimensional lattice with nearest 
neighbor hopping; jumping from site n to site n + 1 with transition probability w^, 
or to site n — 1 with transition probability w~. In this manner, the walker has a 
sojourn probability 1 — — w~ per unit time at site n. The conditional probability, 
P{m,t\n), of finding the walker at site m at time t, given that it was initially at site 
n, satisfies a Markovian master equation for a given realization of the set {wf}- 

We are concerned with the survival and residence probabilities in the finite interval 
D = [—M,L] on the chain. The first is the probability, S'„(t), of remaining in D at 
time t, without exiting, if the walker initially began at site n G D. The second is 
the probability, Rn{t), of finding the particle within the domain D at time i, given 
that it initially was at site n (not necessary in D). The dynamical evolution of both 
quantities follows from the backward master equation [23 

dtFn{t) = wt {Fn+i{t) - F„(<)) + w- (i^„_i(t) - F„(t)) , (2.1) 

where 

F„(t) = J2 Pim,t\ri) . (2.2) 

The survival probability is the solution of the equation l|2.1|) with the initial condition 
Sn{t — 0) = 1, for all n G considering the boundary conditions S_(M+i)i't) = 
SL+i{t) = 0, for all t. On the other hand, the residence probability is the solution of 
the equation (|2.1|) with the initial condition i?„(t = 0) = 1 if n G Z?, or otherwise, 
fulfilling the boundary condition Rn{t) — > 0, for |n| oo for all finite t. 
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Figure 1. Graphical illustration of the definition of the probability transitions. 
The boundary conditions for the first-passage time problem are perfect tramps 
(O ) in the segment extremes. 



Finally, MFPT and MRT can be obtained from the asymptotic limit of the Laplace 
transformed (denoted by hats) survival and residence probabilities, respectively 

Tn = lim Sn{s) , (2.3a) 
T„ = lim i?„(s) . (2.36) 

s— ►O 

2.1. Mean first-passage time 

In the first-passage time problem the random walker only jumps inside a finite interval 
with absorbing ends, as shown in figure ^ The boundary can be simply modeled by 
setting tw^^^^^j^^ ~ '^Z+i ~ (i-^-i the walker cannot jump back into the interval 
once it has been tramped on — (Af -ff)ori + f). From the Laplace transform of the 
evolution equation H2.f|l . using that S'„(0) = 1 for all n G D, and taking the limit of 
equation (|2.3a|l . we obtain the corresponding equation for the MFPT 

u;+ (r„+i-T„)+w- (r„_i-r„) = -f, VneZ?. (2.4) 

This equation must be suplemented with the conditions: r-_(M+i) = ^L+i = 0. 
Equation (|2.4(l is a three-term recursion formula. To get a two-term recursion relation, 
which is simpler to analyze, following Gardiner |14j we make the substitution 

A„ - - T„ . (2.5) 

This yields the equation w,^ A„ — A„_i = —1, and results in = — 7l and 
A_(Af+i) = T^M- Assuming that w+ =^ for all n E D, from these conditions we 
immediately obtain 

A„ = ^ A„_i - ^ , (2.6) 

Wn WA 

^ A„ = -r_M . (2.7) 

nS-D 

Starting from site — Af and applying recursively 1)2.6(1 . we obtain 

" nr 1 ,1-1 n-1 - 

A„^ n ^T^M-^-^ E n (2-8) 
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Note that the last term only runs for n > ~M + 1. Using (|2.8|) in H2.7|l . immediately 
results in 

L k -\ L ^ L fc-l k-l - 
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Finally, writing T„ = T^m+ ^k and substituting in according to equations H2.8|l 

k=-M 

and l|2.9|l . we get 

1+ E n ^ 



7/1+ / L L , fc-1 

k=-Mj=-M 'J / 1 , 1 TT 



1+ E n 



E y.+ + E E n 

k=-M k k=-M+l k i^-M j=i 



k=-Mj = -M 3 



fc-1 fe-1 



,k=-M ^k k=-M+l ^k ,^_M ]=i ^3 



This expression can be additionally recast as 



n— 1 fc 

1+ E n ^+ / 

_ k=-M j = -M 3 

J-n — r 



fc — \ ^ IT 

W^w \k=-M^k k=-M ^k ^=k+ij=k+l ^3 



n it 

E E ^ E n 5: . 

^fc-=-Af '"fc k=-M ^k ^=k+lJ=k+l / 

for n G Z?, and where the sums whose upper limit is n — a run only for n > —M + a. 
The above equation for the MFPT is an exact expression for quenched disorder, given 
that it contains explicitly the full dependence on the basic jump transitions {w^}. 

Let us substitute, in (|2.11|) . = = Wk for all k. This corresponds to a 
symmetrical random walk and we obtain 



^ _n + M + 1 / ^ 1 ^ L-k' 

^—^ 1/17, ^—^ 



L + M + 2y^^^^Wk ^^^^ Wk 

11-2 



E E ^ ^ 



Wk ^ — ' Wfc 

\k=-M k=-M 



2.2. Mean residence time 

In the residence time problem, the walker jumps on the unbounded chain, but we 
compute the probability of finding the walker in the finite region Z), as shown in 
figure |2] Particularly, we are concerned with the mean time that the walker spends in 
D. From the Laplace transform of the evolution equation 1)2. f|l and taking the limit 
of equation l)2.3&|l , we obtain the corresponding equation for the MKT 

w+ (t„+i - r„) + (r„_i - Tn) = -Rn{t = 0) , (2.13) 
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Figure 2. Schematic representation of the residence time problem. The region 
of interest is the segment between the vertical dashed lines. 



where i?„(t = 0) = 1 if n G I?, or otherwise. For this problem, we assume the 
presence of a global bias that points to the right, i.e., w~ /w^ < 1 for all site j. 
Thus, equation H2.13|) must be supplemented with the condition: lim„^_|_oo Tn = 0. 
Equation l|2.13|l is also a three-term recursion formula. Therefore, we make again the 
substitution 

r„ = T„+i - T„ . (2.14) 

This yields the equations 

r„_i = r„ for n < -M , (2.15a) 



w„ 

r„ = ^r„_i forn>L, (2.156) 
'W~ 1 

r„ = ^r„_i -r ioT-M<n<L. (2.15c) 

We are considering residence in a finite region, then we must take the boundary 
conditions: lim„_^±oo r„ = 0. 

Applying recursively equation (|2.15Qf) . we get 

-(M+l) ^+ 

Tk= TT ^r_(M+i) forfc<-(M + 2). (2.16) 
w ■ 

Given that we have assumed Wj^/wJ > 1, the condition limfe__>_oo = 0, imposes 
that r_(Af_|_i) = 0. Therefore, we immediately obtain = for A; < — (Af +1). Thus, 
from equation H2.14|) we get 

Tfe = r_M for k < -M. (2.17) 

In a similar way, from equation l|2.156|l . we get 

rfe= TT forfc>L + l, (2.18) 



and the assumption Wj /wj' < 1 guarantees that limfe^+oorfe = 0, for any arbitrary 
Fi. On the other hand, from equation H2.15c|l . using that F_(jv./+i) 0, we obtain 



F-M = ^ , (2.19a) 

"" -M 



W 
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fe-i fc-i 



fc-1 



^-V^ff^ for -M<fc<L. (2.19&) 



The boundary condition for the MKT when the bias points to the right, 
hm„^+oo Tn = 0, and equation (|2T4|I . allows us to write t„ = — X^fc^n ^k, and using 
equations H2T8|I and (|2T96(I for fc = L, we obtain 



,,.+ 



EinS^ En^^ forn>L, (2.20) 



and 



TL 





E i n 1+ E n ■ i^-^d 



i=-M « j=i+l J / \ k=L+l ]=L + 1 3 



From equation H2.14|l . we can also write t^m — tl ~ X]fc=-Af^fe' using 
equation (j2T96|) . we get, after introducing a change in the order of the sums, 



L-l , L-2 L-l 

T-M ^TL + 



L-l L-2 L-l j - 

E E E n ;;;t- (2-22) 

fc=-A/ fc=-J\/ i=k+l j=k+l 3 

Finally, writting r„ ~ t^m + X]fc=-A/ Tfe, using equation lj2T9b|l . and doing the same 
change in the order of the sums, we obtain 

n— 1 -J n — 2 n— 1 i ^ — 

rn=r-M- — n ^ for-A/<7i<L. (2.23) 

k=-M k=-M j=k+l ^3 

Equations (|2.17|) and 12.2UII - 12.23II arc the exact expressions for the MRT for quenched 
disorder, with the full dependence on the set of transitions {wf}- 

For the unbiased random walk, the residence problem is not defined as can be 
seen from equations H2.20|l - H2.21|) . where the series diverges in the symmetrical case, 
i.e., w'l =w~ . 



3. Homogeneous chain 

The homogeneous biased chain corresponds to the case w'^ — a, wj — b. Introducing 
the bias parameter 7 = b/a, equation (|2.11|l can be easily written as 

L + l-n L + M + 2 7" - 7-^+1 , , 

^-^W^' ail -7) 7-(^^+i)-7^+i for-A/<n<L. (3.1) 
The study of the drift and diffusive regimes of the MFPT follows from equation (|3.1() . 
Additional information about the survival probability in the homogeneous chain was 
reported in references '5? and |24| . 
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For the MKT in the homogeneous chain, we get 
( L + M + 1 



1 

- < 

a 



for n < — Af 

1 -7 
L-n 1 - 

-; + -77^ for -M <n<L (3.2) 

1-7 (1-7) 

1 _ L+M+l 

,2 7"""^ for n > L , 
(1-7)2 

where < 7 < 1. An expression for the residence probabihty in the homogeneous 
chain and the Umit regimes of equation H3.2|) were given in reference |24) . 

4. Weak biased disordered chain 

Two prototype problems are defined in disordered one-dimensional systems. The first 
one, the bond disorder or random-barrier problem, corresponds to the situation in 
which the transfer rates associated with the bond between two neighbor sites are 
symmetric, = In this kind of problem, the statistical properties of the 

system are given by the distribution of the random variable which describes the 
transition rate in each bond. The second class is the site disorder or random-trap 
problem. Here, the hopping probabilities associated with each site are symmetric, 

— w~j^i |21| . and the distribution of hopping rates in each site characterizes the 
model of disorder. Physical realizations corresponding to each class of problems were 
summarized by Alexander et al jSSj. 

In what follows, we consider site disordered chains. We assume that the hopping 
probabilities are strictly positive random variables, chosen independently from 
site to site and identically distributed. However, we are particularly interested in 
the effects of bias in the chain by external fields. Therefore, we admit that the site 
transition probabilities are not necessarily symmetric in the sense that 7^ 'w~. 
Thus, our first random-site model is defined by w'^ = a-l-^j, wj — b + S^j, where a and 
b are positive constants, and {^j} are taken independent but identically distributed 
random variables with {^j) — 0. This form of jump transitions involves an ordered 
biased background with an added random medium. The strength of the bias is given 
by the ratio between a and b and the disorder is characterized by the distribution of 
variables We introduce the parameter e for bias strength by b/a — I — e. This 

selection of parameters allows us to focus our attention in the small bias limit and 
to study the transition to the symmetric diffusive behavior. For practical reasons, 
we can alternatively introduce this additive model by the constrain wJ — — ae, 
where the random variables {w'j'} are taken independent and identically distributed. 
Additionally, wj' > ae for e > 0, otherwise wf is strictly positive for e < 0. Thus, we 



(4.1) 



3 ' J 

can write 

wJ _ ae 

We will find it useful to define the quantities (3k = (^(l/w^)'''^, which we assume finite 

for all fc > 1 . This class of disorder is known as weak disorder |21l 1221 12HI . We 
also define a measure of the fluctuation of the disorder hy J- = (/32 — Pi) / 01- In this 
manner, we can write 

/ -±^\^pl{l^:F5,k). (4.2) 
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In the next step, we will average over the realizations of disorder the corresponding 
expressions for MFPT and MRT, up to first order in the bias parameter e. For this 
model, in both problems we need to use the following expansions 

-A- w~ 1 
wj ^-^ wj 

J— a J J—ct J 

k—a j—Oi 3 k—oc j—a J 

E ^ E n Jiii- E E E E -f;;;^- (4-3^) 

fc=-A/ k k=k+l j=k+l i k=-M k k=-M i=k+l j=k+l k j 

We must note that in the last term of equations (|4.3c|l and Ij4.3djl j ^ k. Therefore, 
the average of these expressions do not involve the fluctuation of disorder. However, 
when we average the equation (|2.11|) . expanding up to first order in e and using 
equations H4.3flf) - (j4.3d|) . we obtain 

(r„) ~ 60(71, L, M) Pi + [61 (n, i, M) + 92 (n, L, M) T\ a0ie, (4.4) 

where the functions Oo(n, 8i(n,i,M) and 82(71, L,M) are defined in the 

appendix. Here, the fluctuation of disorder T is introduced by the expansion of the 
denominator of expression H2.11|l . Following the calculations given in the appendix, 
we obtain 

eo{n,L,M) = ^ '-^ (4.5a) 

@i{n,L,M) L-M + 3-2n 



eo(n,L,M) 6 
62(71, L,M) _ 2L + MH 



eo(n,L,M) 3(L + M + 2) 
Therefore, the averaged MFPT results in 

(L + 1 - n){M + 1 + 71) 



(4.5&) 
(4.5c) 



(Tn) 



, /L-M + 3-271 2L + M + 3-7i_, ^ 
6 + 3(£ + M + 2) -^'"^^^ 



(4.6) 



Taking the limit e ^ in equation (|4.6|) . we arrive to the known expression for the 
MFPT of a homogeneous symmetrical chain with the effective coefficient /3j~^. Note 
that, for this class of disorder, the asymmetry in the hopping transitions links the 
strength of the bias with the fluctuation of disorder. 

On the other hand, for the residence problem we need to evaluate the expression 

\k=n j=L+l J I i=0 ''^ 
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which is vaUd for n > L + 1. Then, taking < e < 1, so that the bias field points to 
the right, and using equations (|4.3e|) . (j4.3d|l and H4.7|l . from equation H2.20|l - H2.23|) we 
obtain the averaged MRT 

{1 for 71 < -M 

(n + M)(n + M + l) ^ , , , 

2iLlM + l) "^^^ for-Af<n<L (4.8) 

l-[n- (L- M)/2] a Pit for n > i . 

This expression is equivalent to the diffusive regime of a homogeneous chain. 
Therefore, it can be obtained from equation (|3.2|) taking 7 = 1 — e and expanding 
up to first order in e. The effect of disorder appears as a renormalization of the 
parameters a and e according to a — > and e ^ a Pit |24j . Note that, in the limit 
e — ^ 0, from equation (|4.8I) we obtain a divergence. This is a direct consequence of the 
fact that the MRT is not a defined quantity for unbiased chains. 

We remark that the results given by equations H4.6I) and H4.8|l are valid for weak 
disorder. When the quantities j3k are not finite, we have to deal with a perturbative 
approach for analyse the behavior of the averaged MFPT and MRT |24) . 

5. Multiplicative model 

Our second asymmetric random-site model is defined by w'^ = a r/j and wj = br/j, 
where a and b are positive constants and {tjj} are taken independent but identically 
distributed positive random variables with (r/j) — 1. Thus, we obtain w~ /wj' = 
b/a = J. It is immediate to recognize that the averaged expressions for MFPT 
and MRT are obtained replacing l/w^ by /3i = {1/ {arjj)) in all the terms of 
equations H2.11|l and (|2.2(J|) - (|2.23|) . Therefore, the averaged formulae for this model of 
disorder corresponds to the homogeneous case, given in section^ with the substitution 
a — > /3f ^ . We must stress that the averaged expressions obtained by this way are exact 
expressions for all value of the bias parameter 7. Particularly, taking 7 = 1 — e, up to 
first order in e, we obtain for the averaged MFPT 



irr. , (L + l-n)(M + l + n) 



L- M + 3-2n 
' + 6 ' 



(5.1) 



2f3i 

As a notable remark, we observe that there is not coupling between the bias e and 
the fiuctuation of disorder J-' for the multiplicative model. We emphasize that the 
multiplicative asymmetric disordered model cannot be easily worked out by using 
FEMA |21l 122 or in general using a perturbative method . 

6. One reflecting boundary 

In this section, wc consider a reflecting boundary condition at the left extreme of the 
interval. The reflecting boundary can be modeled by setting wZj^i — = 0. 

Without lost of generality, we take M = 0, therefore our interval of interest is [0, L]. 
From equation (|2.1H) . taking AI = and Wq = 0, we immediatly obtain for the MFPT 

L -. L — 1 -.Li — 
fc=0 ^fe k=0 ^fe i=k+lj=k+l '^3 

Ti=To-^, (6.1&) 
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E 

fc=0 



1 



T„=T„->;^-^-L ^ n ^ ior2<n<L.{6.1e) 

k=0 i=k+l j=k+l j 

Equation (|6.1a|l was reported for discrete time random walks by Murthy and Kehr \V6\ . 
The main effect of the reflecting boundary is the disappearance of the denominator 
in the expression of the MFPT. For the homogeneous (nondisordered) case, for which 
b and 7 



ri-l 

E 



w: = a, 



b/a, we obtain for < n < L that 
1 



a(l-7) 
The small bias expansion (7 — 



L+l-n- 



7 



7 



L+2 



1-7 

1 — e) of equation Ht).2|l is given by 



T 



(£ + l)(-L + 2) 1) 
2a 



n{n + 1) 
3 ' 3 (L + 2 + n) 



(6.2) 



(6.3) 



In the guideline of sectional we can compute the averaged MFPT for a walker in 
a chain with reflecting-absorbing boundary conditions and additive disorder. In this 
manner, we obtain up to first order in e 



(Tn) 



(£ + l)(£ + 2)-n(n + l) 



n{n + 1) 



3(i 



a Pie 



(6.4) 



Equation (|6.4|l is exactly the same as H6.3|l for weak biased nondisordered chains if 
we set a [3^^ and e — > a Pie. Strikingly, for reflecting-absorbing boundaries the 
fluctuation of disorder is not present in the averaged MFPT. From equation H6.2|l 
and (iniSIl, the same argument presented in section fallows us to reckon the averaged 
MFPT for a disordered chain under multiplicative disorder, substituting a by 

The presence of a reflecting boundary at the left of the interval has not effect on 
the RT problem when the bias points to the right. This fact can be easily seen from 
equations H2.20(l - (|2.23(l . where we found that MKT's expressions do not depend on 
the hopping transitions wZj^.j and w~^^j^.j_^_^y However, if the reflecting boundary is at 
the right of the interval, when the bias points to the right, the MKT diverges as is 



expected. This fact can be seen from equations (|2.21|I - H2.23|I taking 



w 



the limit 



(sec figure I2I1. 



L+l 



and 



7. Concluding remarks 

We have presented an algebraic method for calculating MFPT and MKT of a random 
walk on one-dimensional lattices for quenched disorder. The starting points are the 
one step equations (|2.4() and H2.13() . We have obtained the exact solution (|2.11() for 
the MFPT and the exact expressions ifTTTIl and ((22011- 1223) for the MKT. Also, we 
have considered a reflecting boundary in the sectional Two models for site disorder 
in the chain were considered, namely, additive and multiplicative. The expressions 
for MFPT and MKT were exactly averaged for both kinds of disorder. The main 
difference between these models is in the coupling of the inverse moments of jump 
transitions and the bias parameter in the averaged quantities. For the additive model 
of sectional the terms to first order in the bias parameter are proportional to a /3i e [see 
equations H4.6(l . H4.8|l and H6.4|l ]. Moreover, for absorbing-absorbing boundaries, the 
bias links the strength parameter e with the fluctuation of disorder J- in the averaged 
MFPT [see equation H4.6|l ]. On the other hand, for the multiplicative model, the terms 
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to first order in the bias parameter are proportional to a e and there is not any kind 
of hnking between e and T in the averaged MFPT. 

In this work we have considered models of site disorder. The pair of transition 
probabilities w'^ and w~ , associate to the site j, has the same distribution through the 
random variable or r]j . In addition, the exact equation (|2.10() allow us to consider 
models of bond disorder. In this case, the pair of jump probabilities w'j' and 
associate to the bond between the sites j and j + 1, has the same distribution. For 
tackle this problem, we only need to rewrite the sums in the prcfactor in the following 
manner: 

E n ^-»:m^+ i: ^ n ^ 

Similar arranging can be made with the residence time expressions to deal with bond 
disordered problems. 
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Appendix 

In order to perform the average over disorder for the additive model of section 0] we 
need previously to evaluate the following functions: 

/„(.,M,. g = fe±il)<|±^, (A.1, 

A(.,M). g E E i=(- + ^^)(- + ^^;^)(- + ^'^ + ^), (A.2) 

X — 1 X i 

h{x,M)^ E E E -^.fc^O' (A.3) 

k=-M i=k+l j = k+l 

k^-M i^-M j^-M 
L—1 X i 

Mx,L,M)^ ^ ^ {L-k) 



{L + M) {L + M + 1) (a; + M + I) (x + M + 2) 
4 

{x + M + l){x + M + 2) 



(A.5) 



U{x,L,M)= 2^ 2^ 2^ Sjk = - (A.6) 
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L — l X i 

U{x,L,M)^ E E E (L~k)6,k 

k=-M i=-M j = -M 

_ (3L + 2M - x) {x + M + l){x + M + 2) 

~ 6 ■ 

Using the functions defined above, the functions Qo{n,L,M), Qi{n,L,M) and 
Q2{n, L, M), that appear in equation H4.4|l . are defined by 

n + M + 1 

[(L + M + l) + /o(L,M)] 



(A.7) 



eo(n,i,M) 



L + M + 2 
-[(n + M) + /o(n-l,M)], 



(A.; 



ei{n,L,M) 



n + M +1 
L + M + 2 

1 



1 



{ML,L,M) + UiL,L,M)) 



62(71, L,M) 



n + M +1 
+/i(n-l,M), 

n + M +1 
L + M + 2 

1 



L + M + 2 

iUn - 1, L, M) + U{n - 1, L, M)) - /i(L, M) 



1 



{h[L,L,M) + U[L,L,M)) 



L + M + 2 
ihin-l,L,M) + Mn-l,L,M)) 



(A.9) 



(A.IO) 



n + M + 1 

Equations H4.5afl - H4.5efl follow from equation (|A.8|l - ljA.10p . respectively, replacing the 
functions fi by their explicit expressions. 
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